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Interacting quintessence from a variational approach 
Part I: algebraic couplings 
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We present a new approach to build models of quintessence interacting with dark or baryonic 
matter. We use a variational approach for relativistic fluids to realize an effective description of 
matter fields at the Lagrangian level. The coupling is introduced directly in the action by considering 
a single function mixing the dynamical degrees of freedom of the theory. The resulting gravitational 
field equations are derived by variations with respect to the independent variables. New interesting 
phenomenology can be obtained at both small scales, where new screening mechanisms for scalar 
fields can be realized, and large scales, where one finds an original and rich class of interacting 
quintessence models. The background cosmology of two of these models is studied in detail using 
dynamical system techniques. We find a variety of interesting results: for instance, these models 
contain dark energy dominated late time attractors and scaling solutions, both with early time 
matter dominated epochs and a possible inflationary origin. In general this new approach provides 
the starting point for future in depth studies on new interacting quintessence models. 
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I. INTRODUCTION 

The universe is currently experiencing an epoch of accelerated expansion. During the last fifteen years, since it was 
first suggested by type-la supernovae surveys Qi, this important conclusion has been confirmed by cosmological 
observations of ever increasing precision, which includes measurements of the cosmic microwave background (CMB) 
ii, the Hubble constant Q, baryon acoustic oscillations @ and again type-la supernovae Q- Unfortunately, 
despite all these inputs from astronomical observations, on the theoretical ground we still lack a fully satisfactory 
explanation of this phenomenon. The standard cosmological model, the simplest model capable of fitting all the 
present observational data, assumes the existence of both a non-vanishing cosmological constant and a cold dark 
matter component. The former is introduced to produce the required accelerated expansion at late times, while the 
latter is postulated in order to increase the amount of structure formation needed to be in agreement with observed 
cosmological structures at both large and small scales. The cosmological constant, however, suffers from some profound 
theoretical issues arising from its extremely small measured value if compared with predicted values by (quantum) 
field theoretical considerations; see @ for a recent review. 

In order to solve these issues, or at least to alleviate them, it has been proposed that the current cosmological 
acceleration might be due to the effect of some dynamical field, mimicking the properties of a cosmological constant 
at late times. To agree with current observations such new component must be invisible to electro-magnetic radiation, 
which is the reason why it has been called dark energy, in analogy with dark matter. The fact that both dark matter 
and dark energy cannot be detected through the visible sector could imply new interesting phenomenology. In fact, 
as long as we are unable to probe these dark components but through their gravitational effects, one cannot exclude 
the possibility of non-gravitational interactions between them. Following this reasoning, one is thus naturally led to 
consider cosmological theories where interactions are possible. Note however that if both dark matter and dark energy 
can be measured solely by their gravitational influence, then the physical consequences of their interaction might be 
observationally indistinguishable from similar effects arising in other cosmological models without interaction, as for 
example in modified gravity and warm dark matter models (MU. Nevertheless it is important to study different 
interacting dark energy models not only to search for new interesting phenomenological features, but also to support 
the observational efforts to focus on possible distinctive signatures of such models; a tentative review on these models 
has recently been proposed in [I^ (see also M)- 

The simplest dynamical models of dark energy identify this mysterious component with a minimally coupled canon¬ 
ical scalar field (p. Dark energy theories of this kind are generally known by the name of quintessence and are quite 
popular in the literature because the y y e sufficiently simple to handle and sufficiently complicated to produce non¬ 
trivial dynamics; for reviews see e.g. [131 [l4|. 

On the other hand, the cosmological evolution of (cold) dark matter is usually described by a perfect fluid with 
energy density p and vanishing pressure p = 0. Simple models of interacting dark energy can thus be constructed by 
coupling the scalar field (p to the energy density p of the dark matter fluid. Such a coupling is usually introduced at the 
level of the cosmological field equations adding a term mixing the quintessence and dark matter equations of motion. 
Specifically in the presence of a general coupling between the two dark components, the background cosmological 
equations are commonly expressed as 

(p -f p^) , p -f 3iJp = —Q P(fi -\- 3H (p0 -|- pifP) = Q , (1) 

where H is the Hubble rate, p^ the energy density of quintessence, its pressure and = SttG. We set c = 1. 
An over-dot denotes differentiation with respect to cosmological time and the variable Q quantifies the rate of energy 
exchanged in the dark sector and generally depends on both cp and p (and possibly their time derivatives). If Q > 0 
the energy flows from dark matter to quintessence, while if Q < 0 the energy transfer is in the opposite direction. 
The dependence on p and p of Q is a phenomenological assumption which, together with the form of the scalar field 
potential, determines the interacting quintessence model at hand. Throughout the last fifteen years an impressive 
amount of possible choices has been studied, the most common being Q oc pp, which can easily be motivated in the 
Einstein frame formulation of scalar-tensor theories [isl - fisjl . We cannot discuss in detail this vast amount of literature 
and the interested reader is referred to [T^ . [l^ and in the references therein which provides a general overview of this 
subject (see also [20l - l^ for some recent works). 

In the present paper we will consider a new approach to build models of quintessence interacting with dark matter. 
The coupling will not be added at the level of the cosmological field equations as in Eqs. (Cl), but it will be directly 
introduced into an action characterising both quintessence and dark matter. Variations with respect to the scalar 
field p and the matter degrees of freedom will then provide the coupled equations of motion of dark matter and 
dark energy. We will assume that dark matter can still effectively be characterised by a perfect fluid with vanishing 
pressure even at the Lagrangian level. This implies that a Lagrangian description for perfect fluids must be adopted 
in order not only to obtain the correct expression for the dark matter energy-momentum tensor, but also to define its 
coupling to the quintessence field. 
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There are different Lag rangian approaches to perfect fluids that have been developed during the years; see e.g. Hi- 
[ 2 ^ . Pourtsidou et al |30l| have recently considered a Lagrangian approach to a scalar field coupled to a fluid based on 
the so-called pull-back formalism . In this framework they constructed three new models of interacting dark energy 
with completely new features with respect to the models defined by the standard coupling of Eqs. ©• Subsequently 
they analysed their impact on the cosmological observables showing new interesting phenomenology capable of being 
tested against observational data. The same approach has also been used to study theories of interacting cosmological 
fluids in the context of effective field theories . 

In this paper we will develop the Lagrangian formulation of quintessence coupled to the dark matter fluid. Instead 
of using the pull-back formalism, we will employ the framework with Lagrangian multipliers outlined by Brown in 
[ 2 ^, which is easier to handle and will lead to new phenomenology at both Solar System and cosmological scales. Our 
philosophy will be to fully embrace the variational methods of j28l[ and in particular we will discuss the most general 
actions that one can build out of all the dynamical degrees of freedom associated with the matter fluid. A similar 
approach has already been studied for dark matter non-minimally coupled to gravity where the fluid 4-velocity is 
coupled to the Ricci tensor and the fluid energy density to the Ricci scalar [s^, ■ Here however we will not modify 

the gravitational sector, which will be described by standard general relativity, but we will limit the interaction of 
the matter fluid with the quintessence field. The couplings obtained in this way will first be applied to describe the 
phenomenology of dark energy at small scales, where models presenting environment dependent interactions will be 
able to efficiently screen any scalar field effect and thus to satisfy all Solar System constraints. Then the background 
cosmology obtained within this framework will be studied and the dynamics of different interacting quintessence 
models will be analysed in depth. 

The paper is organised as follows. In Sec.|n]we will briefly review the variational set up for relativistic fluids exposed 
in [ 2 ^. In Sec. Hill this variational approach will be employed to build theories of a perfect fluid interacting with a 
scalar field. The most general algebraic coupling will be considered, deriving the equations of motion and discussing 
the conservation equations. In Sec. HYl this new way of coupling a relativistic fluid with a scalar field will be used to 
build screening mechanisms for dark energy capable of passing every experimental test and Solar System observation. 
Sec.lv] will then be dedicated to cosmology: the cosmological equations will be obtained and the evolution of the 
universe will be analysed with dynamical system techniques. Finally we will draw conclusions and discuss future 
perspectives in Sec. ED 

Notation and conventions: Unless otherwise specified we will assume standard general relativistic notation with 
signature (—, -I-, -I-, -I-) and Greek indices running from 0 to 3. Sometimes the comma notation for partial derivatives 
will be used: for example Units where c = h= 1 will be employed together with = Mp ^ = SttG. 


II. LAGRANGIAN FORMULATION OF RELATIVISTIC FLUIDS 


In this section we outline the variational approach to relativistic flnids following closely the formulation of Brown. We 
will only introdnce the Lagrangian and derive the equations of motion while skipping the detailed discnssions regarding 
the thermodynamical features that can be obtained in this context. The reader interested in more information is 
referred to 

Within Brown’s framework the Lagrangian for the relativistic fluid can be written as 

C-M = —\/—g p(n, s) + , (2) 

where g is the determinant of the metric tensor and p is the energy density of the fluid. We assume p(n, s) to 
be prescribed as a function of n, the particle number density, and s, the entropy density per particle, ip, 9 and I3a 
are all Lagrange multipliers with A taking the values 1, 2, 3 and a a are the Lagrangian coordinates of the fluid. The 
vector-density particle number flux is related to n as 

J>^ = ^gnU\ |T| = ^ 

where is the fluid 4-velocity satisfiyng = — 1. 

The independent dynamical variables which have to be considered in the variation of the Lagrangian @ are ^ 
J^, s, ip, 9, Pa and a^. 

Variation with respect to the metric g^'' gives the following energy-momuntum tensor 


Tfj,u = pU^Uv 



{gfj.Li + U^Uu) , 


( 4 ) 
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which can be rewritten in the standard perfect fluid form 

= pgti^ + {p + p)Uf,U„ , (5) 


once the pressure p has been identified with 



dp 

(6) 

Variation with respect to the other variables gives 



p p^fj. F s9^^ Paoi"^h — 0, 

(7) 

s 

0, 

ds 

(8) 

p 

II 

O 

(9) 

9 


(10) 

/3a 

J^aj, = 0, 

(11) 


iJf^PA),>.=0, 

(12) 

where 

pFp dp 

^ = = A ’ 

n on 

(13) 


is the chemical potential. Eqs. and OT stand for the particle number conservation constraint and the entropy 
exchange constraint, respectively. These can be rewritten also as 


V^(nC/'^) = 0 and V^(s n C/'") = 0 , 


(14) 


where is the covariant derivative with respect to The scalar fields ip and 0 are thus Lagrange multipliers which 
impose these two constraints. In the same manner Pa are three Lagrange multipliers which are needed to restrict the 
fluid’s 4-velocity vector to be directed along the flow lines with constant as specified by Eq. (El)- The meaning of 
Eq. ([5]) can be understood by recalling from the first principle of thermodynamics. If the particle number is conserved 
one can identify the temperature T with 


1 dp 
n dn 


(15) 


Erom Eq. ([8|) we thus have that 9 behaves as a potential for the temperature 


T = U>^9^^. (16) 

Einally Eq. (TT^ determines the dynamics of the Lagrange multipliers Pa, while Eq. o is known as the potential 
representation of the fluid’s 4-velocity. Combining this equation with Eqs. (|H]) and (HU we find 


F = p-Ts = Uf^p^^, (17) 

where F is the chemical free energy. We thus obtain an analogy with Eq. (1161) implying that is a potential for the 
chemical free energy. 

It is possible to show that the fluid field equations 


v^r^" = o, (18) 

are indeed implied by Eqs. dU-dlU. To prove this one can follow the calculations outlined in Appendix assuming 
no interaction with the scalar field; see also [2^. Here we just mention that contracting Eq. dT51) with the projective 
tensor 




(19) 


gives the equations 






p + p 


dX 


( 20 ) 
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which reduces to the usual geodesic equation for a pressure-less fluid (dust); see for example [sj. 

To conclude this section we mention that once equations o-dni) are imposed, one can obtain the on-shell Lagrangian 
by substituting them back into ([2]) . This procedure yields 

>Cm = -V^P = (on-shell), (21) 

with the last equivalence holding up to total derivatives. This means that the on-shell Lagrangian can be represented 
by either negative energy density, or simply by the pres sure. This result does not hold however if the fluid is coupled 
to other fields or non-minimally coupled to gravity |35l-l37l| . 


III. RELATIVISTIC FLUID INTERACTING WITH A SCALAR FIELD: ALGEBRAIC COUPLING 


The standard way of coupling two interacting matter components in general relativity consists in adding a non¬ 
vanishing current Qf^ to the right hand side of the conservation equations. If and are the energy-momentum 
of these two matter components, then their conservation equations can be written as 


= Q. and ) = -Q, 

while the Einstein field equations are still given by 


( 22 ) 


(23) 


This applies also to the case in which one of the two components is the dark matter fluid while the other is a scalar 
field, namely quintessence. 

In what follows, instead of introducing the interaction at the level of the field equations as above, we will employ 
the variational methods developed in Sec. |TT]in order to define a coupling between a scalar field and a perfect fluid at 
the Lagrangian level. 


A. Lagrangian formulation and field equations 


Taking into account the formulation of Sec. im we can now set up our model where a scalar field (j) interacts with a 
matter fluid. The action we will consider is 


S — J d^X (£grav + Gm + + •^int) , 


(24) 


where the matter Lagrangian Cm represents a perfect fluid and is given in (I2|); the gravitational sector £grav is given 
by the standard Einstein-Hilbert Lagrangian 


r _ ^ D 

^grav — 2^2 ’ 

where R is the curvature scalar with respect to the metric the scalar field Lagrangian is given by 


^0 — —y/^ 




(25) 


(26) 


with V a general potential for cj)', finally for the interacting sector we will consider a general algebraic coupling of the 
type 


Ant =-\/^/(n,s,(/)), (27) 

where f(ji, s, (f) is an arbitrary function which will specify the particular model at hand. Note that / cannot depend 
on the Lagrange multipliers, which are needed to impose the constraints (iHD-dni) for the matter fluid and are not 
supposed to mix with the scalar field. In other words, we are assuming that only the dynamical degrees of freedom of 
the fluid, namely n and s, couple with the scalar field. Note that in the Lagrangian approach of Sec. |n]the particle 
number n is not treated as a fundamental variable to use in the variation, but only as a function composed by 
and according to Eqs. (|31). One can thus also consider a more general coupling between the variables J^, s and 
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ij). However if no derivatives are assumed to enter the interacting term, then the function (1271) still represents the 
most general coupling one can build out of the variables J^, s and (f). In the first part of the present work (namely 
the present paper) we will only consider algebraic couplings of the kind (1271) where no derivatives appear. Interacting 
terms with derivatives will be discussed in the second part [s^. 

Action (|24)) has to be varied with respect to the fields , </), s, ip, 9, Pa and a^. However the variation with 
respect to the last four among these will give again Eqs. (|^- (IT^ and thus will not be repeated. 

We start considering the variation in s which produces the following equation 




1 dp I df 

_L.-I_ L = T + Ti 

ndn ndn 


(28) 


The scalar field 9 represents now a potential not only for the fluid temperature T but also for the “interacting 
temperature” Tint. Variation with respect to yields the equation 

{p + Pint) Ufj. + ^,(1 + + PAOpfj = 0 , (29) 

where pint = df /dn has been defined in analogy with the chemical potential (1131) . This equation modifies the potential 
representation of the fluid 4-velocity d?]). Using Eqs. (fTTl) . (1^ and (0^1) we can show that 


=F + Tint, (30) 

where Tint = Mint — sTint- This equation generalizes Eq. (II3 and tells us that ip is now a potential for the chemical 
free energy of the fluid plus its interacting counterpart. In general Eqs. dm)-®, together with Eqs. (©-dni), define 
how the thermodynamical properties of the fluid change due to the interaction with the scalar field. We will not 
discuss these relations further but focus on the spacetime dynamics only. 

At this point we perform the variation with respect to the metric tensor which results in the following Einstein field 
equations 


n — (rp I rp{4,) . rT-.(int)\ 

— rv pi/ j J 

where is the fluid energy-momentum tensor given in dS]) and 




dpcj)d„4> - ppi, 


^df,cl>d>^cl> + V{(j)) 


(31) 


(32) 


is the energy-momentum tensor of the scalar field </>. The interacting energy-momentum tensor can be written 

as 


= Pint Ppu + {Pint + Pint) UpU„ , 


(33) 


where we have defined 


Pint = f{n, s, (j)) 


and 


Pint = n 


df{n,s,(l)) 

dn 


f{n,s,(l)), 


(34) 


as the interacting energy density and pressure. Finally the variation with respect to the scalar field yields the following 
modified Klein-Gordon equation 


ap- 


d(j) 



(35) 


where □ = V^V^. 


B. Conservation equations 

It is interesting to reformulate the field equations above as in the standard approach to coupled matter component 
in general relativity, i.e. as defined by Eqs. (I22|) - (I23I) . For this purpose one can dehne a new energy-momentum tensor 
for the fluid as T^jy = T^,y -I- T^™*^ with the energy density p = p + pint and pressure p = p + Pmt- In this case the 
Einstein field equations (|^ become 

Gpi/ = ^ (r^. + T^«) , 


( 36 ) 
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resembling Eqs. (l23l) . Using the Klein-Gordon equation (1351) and the energy-momentum tensor (l3^ . the conservation 
equation for the scalar field can be written as 


= Q, , (37) 

which shows that the scalar field is not conserved due to the interaction with the fluid. 

To prove that also the fluid energy-momentum is not conserved, it is better to split its conservation equation into 
the parallel and perpendicular components to the fluid flow: 

, (38) 

where is defined in Eq. dni). Then using Eqs. (I^- (IT^ and Eq. (|29)) one can show that (see Appendix El) 

and . (39) 

Inserting this back into Eq. (1551) gives 


= (40) 

This shows that the Lagrangian formulation we are considering in the present paper can be mapped back into the 
standard relativistic approach defined by Eqs. (I22l) - (l23l) . Note however that the energy density and pressure of the 
fluid differ from the uncoupled case in as much as they now depend also on the scalar held through the interacting 
term. In particular even if p = 0, e.g. in cold dark matter applications, the pressure p is not expected to vanish. As 
we will see these new features will give rise to a new interesting phenomenology at both small and large scales. 

The argument above can also be repeated directly with which does not depend on the scalar field. This provides 

I/.VaT^" = 0 and = -2nt/^V[A (t/^]Pint) , (41) 

where square brackets between two indices denote antisymmetrisation. Eq. (1411) implies that in general 

^ 0 , (42) 

meaning that the uncoupled part of the matter fluid is no more conserved in the presence of the coupling to the 
scalar field. Note however that the longitudinal part of the conserved equation, the one aligned with [/^, is always 
conserved. This will be relevant for cosmology where the isotropy and homogeneity of the spacetime will imply the 
conservations of at the background level. 


IV. SMALL SCALE PHENOMENOLOGY 

Before employing the results developed so far to construct new models of quintessence coupled to dark matter, we 
study their possible applications to Solar System and galaxy scale phenomenology. For this purpose in this section we 
will use the Lagrangian approach developed in Sec. IIII Al to investigate new couplings between the scalar field 4> and 
baryonic matter. In what follows thus p and p will represent the energy density and pressure of matter constituted 
by Standard Model particles. 


A. Fifth force 


We start recalling that because of the Bianchi identity the relation V = 0 holds. We then expect the right 
hand side of Eq. (1311) to be covariantly conserved, giving in this way the equations of motion of the matter fluid. 
Taking the covariant derivative of Eq. contracting with hfn, as given in Eq. m, and using the scalar field 

equation (l35l) . eventually yields 


dUi^ 


+ = 


1 

P+P + Pint + Pint 


P,v + Pint,u 





dX 


(43) 
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which generalises Eq. dsni). For a pressure-less fluid, or equivalently for poirrt particles, this expression reduces to 


dX "" 


1 

P Pint ~t“ Pint 




h''^^ = Z'". 


(44) 


Comparing Eqs. and (EUl) one realises that the interaction between the scalar field and the matter fluid gives rise 
to a fifth force /^. Thus, as one could expect, the motion of the matter fluid is non-geodesic due to the presence of 
the extra force which in general depends on both n and (j). Note however that if does not depend on either cj) 
or n then it vanishes completely, reflecting the fact that there is no effective interaction in such cases. Furthermore, 
according to the definition of 4-force in relativistic theories, is always orthogonal to the fluid 4-velocity: = 0. 

The appearance of this fifth force can give rise to interesting phenomenology at both Solar System and galactic 
scales. For example the non-geodesic motion of luminous matter can be used as an alternative to dark matter in 
explaining the anomalies observed in the galaxy rotation curves. The situation here is similar to the one arising from 
non-minimally coupled theories [s^ , where an extra force depending on the local curvature produces effects similar to 
the ones predicted by MOND theories [4lj. In MOND theories however the modification of Newtonian mechanics at 
galactic distances depends on an effective acceleration qm (or equivalently on an effective length Im) which is usually 
assumed to be constant and thus difficult to adapt to different galaxy profiles. In non-minimally coupled theories 
such effective acceleration depends on the small value of the local gravitational curvature as well as on the matter 
energy density of the environment. This might help in fitting the rotational velocity curves of several galaxies within 
a single model where om varies from one galaxy to the other. Exactly the same argument applies to the scalar field 
interacting model considered here. From the fifth force it is in fact possible to define an effective acceleration 
Qrf, which modifies the Newtonian equations at galactic distances in exactly the same way as om does. The difference 
resides in the fact that depends on the scalar field 0 rather than on the local curvature and thus might predict 
quite different results even for galaxies having roughly the same matter density distribution. 

At Solar System scales the fifth force (IT41) is highly constrained by the experiments which do not show any violation 
on the inverse-square Newtonian law up to parts in 10® [d^. The effects of the interaction between the scalar field and 
the matter fluid must then be negligible at those distances. Fortunately, as we are now going to see, the dependence 
of the fifth force on the fluid dynamical degrees of freedom allows us to build efficient screening of possible deviations 
from the geodesic motion. 


B. Chameleon mechanism 


From the Klein-Gordon equation (l35l) one can realise that the scalar field (j) “feels” an effective potential given by 

Ves = ¥{(/)) +f{n,s,(l)), (45) 

which can be used to screen the interaction between the scalar field and matter at Solar System scales. The effective 
potential (H51) represents in fact a generalisation of the so called chameleon theories where the scalar field couples to 
the matter energy density and acquire an effective mass characterised by the surrounding environment . To show 
this we can consider an interacting term defined by 

f{n,s,cj}) = p{n,s)j{(j)), (46) 


where p is the matter energy density and 7 is a function of (j). With this choice the interacting energy density and 
pressure m read 


Pint = p{n, s) j{4>) and pint = p{n, s) j{4>). 
Next we assume the scalar field potential to be of the “runaway kind”, i.e. to satisfy 


lim V = 0 , 

(p—^OO 


lim y = 00 , 

0 -tO 


dii)v 


lim 

(p—^CO 

lim „ , 


/ 




= 0 


ab-i)y 

i9Z)(»-i) 


for i = 1,2, 3,..., 
for i = 1,2, 3,.... 


(47) 


(48) 

(49) 


Such constraints are usually imposed on the scalar field potential in order to obtain a late time cosmic speed-up such 
as, for example, in the well-known inverse power law potential [d^ 


vw 




( 50 ) 
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where M is a constant with mass units and a. a positive parameter. If we now set 

7 ((?i) = 706 ”'^'^ , (51) 

with P and 70 positive constants, we obtain the effective potential 

Kff = V(</>)+7ope"'5^ (52) 

which coincides with the single matter chameleon potential of [43 | . Of course if we considered more than one matter 
fluid we would have ended up with exactly the multi-component equation considered in (d^ . however for our scopes 
a discussion with one matter fluid is enough. The fact that the Ves has a minimum for Solar System values of p 
implies that at those distances the scalar field acquires an effective mass. If this mass is sufficiently high, then all 
possible interactions with particles of the Standard Model are efficiently suppressed. At cosmological scales though, 
where the matter energy density is extremely low, the effective potential (I52|) reduces simply to V((f>) which, given 
the assumptions (|4^ and (|4^ . can thus drive the late time acceleration of the universe; see for more details. 

From these considerations we have seen that the formulation developed in Sec. lIII Al with the assumption f{n, s, (/>) = 
p{n, s) 7 ((()), can be used to derive chameleon field theories from a Lagrangian approach. Moreover if other interacting 
terms /(n, s, 4>) are considered, the same formulation can also be adopted to build more general screening mechanisms 
for the scalar field. In fact as long as the effective potential (1451) posses a minimum for Solar System values of n and 
s then the scalar field will always appear with a non-zero mass at those scales, which, if large enough, will allow the 
scalar field to pass all Solar System tests. Of course, in order to properly check its viability, any such model should 
undergo an analysis similar to the one performed for the original chameleon theories of . A deeper analysis of the 
possible screening mechanisms that can be obtained with different choices of the interacting term /(n, s, p) will be 
left for future studies. Here, in order to present some of the potentialities that derive from our approach, we will limit 
the discussion on few considerations regarding a particular model. 

First we will assume the function /(n, s, p) to depend on n and s only through p(n, s). This reduces the space of 
possible new models, but allows for an immediate physical interpretation and the analysis can be better compared 
with the standard chameleon paradigm as given by Eq. (Id 6 l) . In any case the reader should keep in mind that from 
the Lagrangian formulation developed in Sec. IIII Al it is possible to construct screening mechanism where the mass 
of the scalar field depends on the local number density n and entropy s. Note that though the effects of n might be 
similar to the one obtained with p, the possibility of creating “entropic screenings” using different values of s is an 
interesting feature arising from our approach which might be worthy investigate in future workt 0 . 

For the particular model we will consider here, we will take a scalar field with vanishing self-interaction, i.e. with 
V{(j)) = 0. Without a coupling to the matter sector such massless scalar field could never accelerate the late time 
expansion of the universe and thus it would never constitute a viable dark energy model. We will however assume an 
interacting term of the general form 


I4ff 


f{p,(j)) = A cosh 




(53) 


where /3s and a are two dimensionless parameters, A is the cosmological constant and Mp = 1 /k = I/VSttG is the 
reduced Planck mass. For any fixed value of the matter energy density p, the mass of the scalar field is given by 


/I ^ 



(54) 


To drive the late time accelerated expansion of the universe, the inverse of this mass, i.e. the Compton wavelength 
A^, should be of the order of the Hubble length 1/Hq at cosmological scales, where the energy density is roughly 
Pc — 10“^®g/cm^. In this manner the effective potential (1531) reduces to nothing but the standard cosmological 
constant for sufficiently large distances and thus it provides an accelerated expansion with the simplest possible 
mechanism. On the other hand if p is sufficiently high around laboratories on Earth, then the interaction mediated 
by the scalar field will be efficiently suppressed because of its large mass. This situation is qualitatively depicted in 
Fig. [T] where it is shown how the potential (1551) changes for different values of p. 

Imposing now the condition = I/LIq Pc one finds 


= H, 


Mp 

VX 


~ j^QlO-|-122a 


(55) 


^ One might also think to relax the condition on the conservation on the entropy density (llOH to probe possible non-adiabatic effects. 
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FIG. 1. Qualitative behaviour of the interacting potential (I53II for different values of p. On Earth (ps) the scalar field is highly 
massive, while for galactic densities (pgai) the mass is lower. At cosmological scales (pc) the mass is so low that the potential 
effectively becomes a cosmological constant. 


which given any positive value of a represents an anomalously large number for the coupling constant /3s. Nevertheless, 
leaving this problem aside, one can now look at the values the mass takes on the Earth and within a galaxy. Using 
the order of magnitude values pE — lOg/cm^, patm — 10“^g/cm^ and pgai — lO^^'^g/cm^ respectively for the energy 
density inside the Earth, in our atmosphere and in the galaxy, one finds 


As an example we consider the values a = 1 and a = 2 which yield 

~ 10“"^ m, A^*™ ~ 1 m, ~ 10^^ m ~ 100 kpc, 

Af~ 10-34 m, A^*“ ~ 10-26 m, A®*^^ ~ lO^® m ~ 1 pc 


(56) 


(a = 1) (57) 

(a = 2). (58) 


In both cases the interactions mediated by cj) are short-range on the Earth and long-range on galactic scales. In 
the case a = 2 the mass on the Earth is so high that any effect of the scalar field is practically undetectable. 
Nevertheless in both cases at galactic scales the scalar field characterises a long-range force which might give rise to 
the interesting phenomenology related to the anomalous rotations of galaxies discussed in the previous subsection. 
The same fifth-force could also produce testable effects at the level of clustering scales which are easy to analyse 
within linear cosmological perturbations theory. From these considerations we obtain that within a single scalar held 
model dehned by the interaction (1551) one could in principle explain dark energy, dark matter and automatically pass 
all Solar System constraints. 

These brief calculations show the potential of the approach developed in this paper as a powerful tool to construct 
efficient screening mechanisms which can also provide desirable effects at both galactic and cosmological scales. We 
will now turn to the study of cosmological applications of the framework of Sec. IIII Al and leave further analyses at 
Solar System and galactic distances, as well as at the linear perturbation level, for future investigations. 


V. COSMOLOGY 


In this section we will consider the evolution of the Universe as a whole using models based on the interacting 
formulation developed in Sec. IIII Al and deriving suitable equations for their dynamics. The fluid matter sector will 
now describe dark matter and thus the interacting term between the scalar field and the fluid variables must be 
considered as an effective coupling between dark energy and dark matter. Note that due to the results of Sec. HVl 
we could also consider a viable coupling to luminous matter if an efficient screening acts at Solar System distances. 
Nevertheless we will only analyse a coupling to dark matter since this dark component dominates over baryonic matter 
in the cosmological matter sector. 












11 


A. Cosmological equations 


We start considering a Friedmann-Robertson-Walker (FRW) line element as required by the cosmological principle. 
The metric will then be determined by 

+ , (59) 

where a{t) is the cosmological scale factor, k = —1,0,1 according to the spatial openness, flatness or closeness and 
is the two-dimensional line element of a sphere. We will also assume that all the dynamical quantities are 
homogeneous, i.e. they depend only on t. In particular we will have that 0, p, n, s will be functions of t only. 
Moreover taking into account coomoving coordinates the perfect fluid 4-velocity becomes simply = (—1,0,0,0). 

The cosmological dynamics is determined by Eqs. (fTdll . (1^ and (l35|l . However in a cosmological setting we will 
always have that Eqs. d yields 

h + 3Hn = 0 and s = 0, (60) 

where the over-dot denotes differentiation with respect to t and H = d/a. These equations tells us that the entropy 

density per particle is conserved through the Universe evolution, while the particle density decays according to 

n oc a~^ , (61) 

which is expected from geometric considerations. The only two dynamical quantities that determine the evolution 
of the Universe are thus a and (j) whose behavior is controlled by Eqs. m and (1551) . In a FRW cosmology those 
equations yield three independent evolution equations, namely the two Friedmann equations 

^ + V + Pint^ , (62) 

- u + Pint) , (63) 

and the scalar field equation 

+ + (64) 

Since at the cosmological scales all the dynamical helds are effectively homogeneous, the right hand side of Eq. (PI) 
vanishes. More explicitly for every dynamical quantity x have x,^, = ~xUfj, and thus = 0. It follows that 

the matter equation of motion are not modified at large cosmological scales 

p + 3H{p + p) = 0, (65) 


3— -h 3H‘^ = 


1 1 ‘ 

'’+U' 


+ 2H + 3H^ = -K^ [p 


ds^ = —dt^ + a(t)^ 


dr'^ 


1 — kr^ 


and, assuming p = wp, the energy density decays exactly as in a non interacting universe 

p oc . ( 66 ) 

This is in agreement with the conservations equations m since in a homogeneous and isotropic cosmological frame¬ 
work the transverse part of those equations vanishes; see below for a derivation using the cosmological equations. 

At late times a 1 the cosmological dynamics will be completely determined by the scalar field (f) since both n 
and p decay and pint and pint will depend on (j) only. If in our model we have that pi„t = 0 at late times, then pint just 
describes further dust (dark matter) and, provided p^ ^ pint at late times, the universe evolution will be determined 
by the form of the potential V{(j)). Choosing a suitable quintessence potential will thus gives rise to a late time cosmic 
speed-up as required by observations. 


B. The conservation eqnation 


As it is not trivial to see that the standard matter conservation equation indeed holds, we will show this explicitly 
in the following. We consider the derivative of (15^ and the combination (15^ minus (155)) . Respectively, they give 


-3h(2^-2h'^ =K^{p + , 


2— -2H = K^(p-|-p 


’ + y'4> + Pint) 

\- (j)^ + Pint + Pint)- 


(67) 

( 68 ) 
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Therefore we get 


P + ^(j) + V'(j) + Pint — —3H(p + p + + Pint + Pint)- 

A few algebraic manipulations give the equation 

P + 3-ff {p + p) + (j) (j) + 3Hlj) -VV' + Pint + 3ij (pint + Pint) = 0. 
Using (El, the last terms involving pint and pint can be expanded as 

^ I ozj/ I „ ^ ^ I ^Pint . I d Pint I , 

-Pint + 3il(pint + Pint) - + -^s + 

dP\nt ! . I o tJ ^ I ^Pint . , Spint 

The first and second term are zero due to (16011 . Now we can rewrite (j70l) in the following way 

dpint ■ 


p + 3iP (p + p) + 9^ U 5H<p + V' + 


dcj) - 


= 0 . 


(69) 

(70) 

(71) 

(72) 

(73) 


The term proportional to (j) is zero because of the equation of motion of the scalar field (1641) . Therefore, we have 
proved that 


p + 3i7 (p + p) = 0 , 


(74) 


in agreement with Eqs. (16511 . 


C. Relation with standard interacting models 


In the following we will clarify the relationship between the variational approach used here and the frequently used 
approach of introducing couplings at the level of the conservation equations. Using the notation p = p + pi„t and 
p = p + Pint one can rewrite Eq. (1701) in the form 


p + 3H {p + p) + <j) (ji + 3H^ + V' 


= 0 , 


(75) 


which is the sum of two conservation equations, that of the fluid with p, p, and the conservation equation of the scalar 
field. Let us now recall the equation of motion of the scalar field (I64|) which we write 

+ = ^ with Q = (76) 

0(p (f) d(p 

so that it takes the standard form with coupling Q, see also Eqs. O- Consequently, Eq. (ESI implies a coupled 
equation for the matter 


^ + 3iJ(p + p) = -Q. (77) 

Note that these equations can equivalently be derived from the general covariant equations dSZI and dsni- 

We are now in a position to link the two different approaches to coupled cosmological models. Our approach is 
equivalent to standard coupled models provided we identify p and p with the usual matter variables p and p in the 
standard approach. In other words, the standard approach uses an effective description of the matter sources present. 
The quantity p, for instance, contains a scalar field dependence and therefore the standard approach does not allow 
the consistent separation into matter and interaction parts without additional ad-hoc assumptions. This is not an 
issue solely concerned within scalar filed models, but it equally arises in any dark energy model interacting with dark 
matter, as recently pointed out in 

At least for scalar fields, we can conclude that every standard coupled model can be derived from the Lagrangian 
formulation presented here. However, it is also clear that the correct identification of matter is required to make this 
approach consistent. 
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D. Dynamical system techniques 

The subsequent analysis of the background dynamics relies on the use of dynamical systems techniques. In order 
to be self-contained, we briefly introduce the most important concepts tailored to three-dimensional systems. A 
three-dimensional autonomous system of differential equations, also called a dynamical system, is given by the three 
equations 


x' = fi{x,y,z), y' = f 2 {x,y,z), z' = f'i{x,y, z), 


(78) 


where the prime denotes differentiation with respect to a suitable time parameter. The three functions fi{x,y^z) do 
not explicitly depend on the time parameter, they only depend on the dynamical variables x, y, z. In technical terms 
this means we are dealing with autonomous systems only. 

A fixed point, critical point or stationary point of the system (17811 is a point with coordinates {xQ,yo, Zq) such that 
fi{xo, yo, Zq) = 0 for i = 1,2,3. This corresponds to a stationary point of a mechanical system where the potential has 
an extremal point and a fictitious particle would remain at rest. In analogy, at the point (xq, yo, zq) the system is at 
rest since all evolution equations are identically satisfied; in principle the system could remain in this state indefinitely. 

Linear stability theory clarifies whether a stationary point is stable or unstable with respect to small perturbations 
away from that point. The basic idea behind this approach is to Taylor expand the three functions fi{x, y, z) around 
the fixed point (xq, yo, zq). Since the functions vanish at the fixed point, the first non-trivial terms in the Taylor series 
will involve the first partial derivatives of the functions. Let us denote Xj = {x,y,z) for j = 1,2,3. It is therefore 
natural to define the matrix of first derivatives 



hj — 1 ) 2 , 3 , 


(79) 


which is the familiar Jacobian matrix of vector calculus (not to be confused with the current J^). The information 
on stability is contained in the eigenvalues of this matrix J evaluated at the critical point {xo,yo, zq). Since J is a 
3x3 matrix, it will have three, not necessarily distinct, eigenvalues. If all eigenvalues of J are positive, we speak of 
an unstable point since all perturbations would grow exponentially. On the other hand, if all eigenvalues are negative 
this point would be regarded stable. If at some fixed point the matrix J contains positive and negative eigenvalues, 
then one speaks of a saddle point. 

One could also encounter a pair of complex conjugate eigenvalues and one additional real eigenvalue. In this case, 
stability and instability will depend on the signs of the real parts of the eigenvalues. Linear stability theory breaks 
down when at least one of the eigenvalues has a zero real part. In this case one has to use techniques beyond linear 
stability theory to understand the dynamics of the system near this fixed point, see for instance |45| . However, these 
techniques will not be required in the following as linear stability theory turns out to be sufficient. 


E. Cosmological dynamics 


In this section we will study the dynamics of a universe described by Eqs. (EUl-dMl) employing dynamical systems 
methods to determine its complete dynamical properties. Similar analyses for dark energy models coupled directly 
through the cosmological equations, i.e. considering Eqs. © for different phenomenological choices of Q, have already 
been performed in several works; see e.g. [^1^1^. We introduce the following dimensionless variables 


2 i^'^P 

G = 


3i72 


x2 = 


6i72 


V^ = 


3772 > 


3772 


(80) 


which straightforwardly generalise the normalised variables usually employed to analyse the dynamics of quintessence 
[ 4 ^. In general, since pint is a function of n and (j) {s is constant because of dini) and will not be considered in what 
follows) which in turn can be seen as functions of p and V respectively, we can always write it as a combination of 
the variables a and y. However, unless the function pint is specifically choosen not to do so, this will require the 
introduction of another variable such as z = 77/(77o -I- 77) with Ho a constant. To directly consider a fourth variable 
z as defined in (1501) simplifies the analysis since it allows us to reduce the Eriedmann equation (1001) to the constraint 

1 = +y‘^ + z , (81) 


which will permit us to replace cr in terms of the other variables. Moreover we consider only y > 0 since we can 
assume V > 0. Alternatively we will see that the system will be invariant under the changes y 1 —> —y, meaning that 
its dynamics will be specular over the xy planes. 
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In what follows we will assume that p = wp with w a constant, called the matter equation of state parameter, 
physically constrained to be between 0 and 1/3. For (dark) matter one has w = 0, however we will keep w as a general 
constant for the moment. We will also assume the potential V to have an exponential form of the kind 


= Fo , 


(82) 


with Vq a constant and A a dimensionless parameter. The exponential form for V is the only one which consents to 
close the autonomous system of equations without introducing another variable. This happens because V' (x V and 
thus derivatives of the potential can be related to y. For different potentials instead a further variable must be added 
in order to take into account the derivatives of the scalar field potential 

With these assumptions the cosmological field equations (l62l) - (l64l) can be rewritten as 


x' = —— (^3x ((in + l)y^ + wz — in + l) + 3(in — l)a;^ — \/6Xy^^ + xA — B (83) 

y' = ~ 2 ^ ~ 1)®^ + 3 ((in + l)y‘^ + inz — in — l) + V^Xx'j + yA (84) 

z' = 2A{z — 1) + 2Bx — 3z ((in — l)x^ + (in + l)y‘^ + w{z — 1)) . (85) 


where a prime denote differentiation with respect to dp = Hdt and we have defined 

K dpint 


A = 


K^Pint 


21/2 


and B = 


^/QH^ d(j) 


( 86 ) 


In order to close the system we must specify the function pint and then compute A and B using (IMl) . If pint is chosen 
accordingly the quantities A and B become functions of x, y, z and the system results to be close. If instead A and B 
cannot be written as functions of x, y, z then another variable has to be added in the dynamics and the dimensions 
of the system increase. 

The following choices do not increase the dimensions of the system further: where a, P and 7 are dimensionless 



Pint 

pint 

A 

B 

Model I 

7 p“ exp(—/I k^) 

[a(m -1- 1) - l]pint 

|[a(ui - 1 - 1 ) — l]z 


Model II 

JKipp 

mpint 

|mz 



parameters. These choices are the simplest ones found by the authors which do not require a dynamical systems 
of more than three dimensions, i.e. for which A and B can be written as functions of x, y, z. Mathematically 
they are simple to analyze and physically they are sufficiently complex to allow us to investigate the new and rich 
phenomenology of the scalar-fluid coupling considered in this work. 

These two models are also characterised by a linear equation of state pint = iHintPint with Wint being constant and thus 
they represent the simplest three-dimensional models one can think of. However, by considering other forms of pint, 
one can construct models where the equation of state iCint itself becomes dependent on the variables x, y, z. One such 
choice is pint = exp(p/po) where po is some constant characteristic density. This gives pint = [(1 + w)p/po — l]pint 
which means icint = (1 -I- w)p/po — 1. Recall that p is related to which satisfies (IM|l and thus introduces the 
dependence on x,y,z. 

Unfortunately, such a model would increase the dimensions of the corresponding dynamical system since mint 
cannot be written solely in terms of the variables x, y, z. A coupling where this can be achieved corresponds for 
example to the rather complicated choice pint = 7 U/ln(p/po) with U((/)) the potential (15^ . In this case one obtains 
mint = —I — (m -|- I)z/( 7 y 2 ) and both A and B can be written as functions of x, y, z without introducing other 
dynamical variables. 

In what follows we will consider only the two simplest choices corresponding to a linear interacting equation of 
state and leave analyses of the more complicated models for future studies. We will label them as Model I and II as 
outlined by the table above. 


1. Model I 

In this subsection we analyse Model I of the dynamical system (1551) - (1551) . Now the two Friedmann equations (1551) 
and (1631) can be written to give the acceleration equation 

^ = ^[-(1 + w) + {w- l)x^ + (1 + m)y^ - (1 + w){a - l)z] 


(87) 
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which can be solved for a at any fixed point (x*, y*, z^,) of the phase space to give 

a oc (t — (88) 

where tp is a constant of integration. This represents a power law solution, where the scale factor a evolves as a power 
of cosmological time t. If the denominator of (1551) vanishes, this forces H to be constant, which corresponds to the 
universe undergoing a de Sitter expansion. If we define the effective energy density and pressure of the system to be 

PeS = P + 2 *^^ + V + Pint (89) 

PeS=P + - V + Pint (90) 

we can then define an effective equation of state paramter Weff, which can be expressed in terms of the dimensionless 
variables (1501) 


PeS 

WeS = - 

PeS 

= w — {w — l)x^ — (1 + w)y^ + (1 + w){a — l)z. 

(91) 

These definitions allow us to write (1551) 

in the simple form 




(92) 

The energy density and pressure of the 

scalar field are given by 



P<P = 2*^^ ^ 

(93) 


p<i> = 

(94) 


and we can also define the equation of state of the scalar field at any point in our phase space 


Pcj, 

W4, = — = 

p<j> + y‘‘ 

Now that we have specified the model, we can find the critical/fixed points of the autonomous system 
These are defined to be the points (x, p, z) of phase space which satisfy 

x' = 0, y' = 0, z' = 0 


(95) 
-(ES]). 

(96) 


If the system is at one of these points, there is no dynamical evolution and the universe evolves according to (1051) . In 
order to exist, the critical points must lie in the three dimensional phase space, which are all points (x, y, z) such that 

x^+2/^+z<l. (97) 


This follows from the assumption cr^ > 0, i.e. p > 0, and the Friedmann constraint (15T1) . Stability of critical points 
is determined by linearising the autonomous system around the critical point under consideration, which leads to 
analysing the eigenvalues of the Jacobian matrix 


Aiij 


dfi 

dxi 


(98) 


evaluated at the critical point, where we have compactly written the dynamical system (1551) - (155|) as x' = fx(x,y,z), 
y' = fy{x, y, z) and z' = fz{x, y, z). If all three eigenvalues of A4 have negative real part the critical point is stable, 
if all three have positive real part the point is unstable, and if some eigenvalues have different signs the point is a 
saddle. 

There are up to eight critical points of the system, depending on the values of the parameters a, /3, w and A. They 
are outlined in Tab. HI 

We have checked explicitly that there are no critical points at infinity. We did this by using the compact variable 
C = arctan(2:) and it turns out that there are no critical points where = ±7r/2. 
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Point 

X 

y 

z 

0 

0 

0 

0 

A± 

±1 

0 

0 

B 

/?(i+™) 

V 2 A 

/i" Y (i+u')(i-'u') 

V 2 A 

0 

C 

A 

a/6 


0 

D 

n p 

Y 3 a(w+l)-2 

0 

1 2/32 

Z{ol-\-oiw—2)'^ 

E 

fz (l + i(j)(l —O') 

0 

3 (1 — o')(l+L(j)(l—ly) 

V 2 P 

2 'P 

F 

fz (1+l(;)q: 

y/&{w+l)a-Z{w+l)'2a'2+2P{P-X) 

A(A—/?)— 3(l+iiu)Q: 

V 2 A -/3 

y/2\\-P\ 



TABLE I. Critical points of Model I. 


Point 

X 

y 

z 

0 

0 

0 

0 

A± 

±1 

0 

0 

B 

/T 1 

V 2 A 

fz 1 

V 2 A 

0 

C 

A 

a/6 

A-if 

0 

D 


0 

1 2/32 

^ 3 

F 

/T 1 

V 2 A -/3 

y/ 3 + 2 P(P-X) 

V 2 \X-P\ 

A(A -^)-3 

(/ 3 -A )2 


TABLE 11. Critical points of Model I with w = 0 and a = 1. 


The case a = 1. 

The dynamics of the system depends on the four constants a, /3, A and w. In order to simplify the analysis we will 
first consider the natural choice a = 1. This particular choice makes the critical point E coincide with the origin O. 
We will also only analyse the dynamics in a matter dominated universe with w = 0, the dynamics of the system for 
other values of w are qualitatively the same and not relevant for models interacting with (cold) dark matter. We are 
thus left with the critical points listed in Tab. El 

Properties of these critical points, including existence and stability can be found in table IIIll There are potentially 
up to seven critical points depending on the values of A and /3: 

• Point O. The origin of the phase space exists for all values of A and P and corresponds to a matter dominated 
universe. In this case Weff = 0 (or more generally WeS = w), so there is no acceleration at this point. This point 
is always a saddle point. 

• Point A±. These two points are dominated by the scalar field kinetic energy, with the effective equation of state 
reducing to that of a stiff fluid WeS = 1- Thus no acceleration is present at this point. These points are either 
unstable or saddle points depending on whether the absolute values of A and /3 are less than VG and 3/-\/6 
respectively. 

• Point B. This point exists for all values of A and p. It corresponds to the usual scaling solution where the 
effective equation of state matches the matter equation of state, yet the scalar field energy density does not 
vanish. Hence the universe behaves as if it were completely matter dominated, yet the energy density of both 
the matter field and scalar field do not vanish. This point is only stable when /3 > 0. 

• Point C. This point corresponds to a universe completely dominated by a scalar field. It exists only for A^ < 6. 
It is stable when A^ < 3 and A/3 > A^ — 3, and a saddle node otherwise. This point is the usual cosmological 
accelerated expansion driven by a sufficiently flat scalar potential: A^ < 2. 

• Point D. This point exists for all values of A and /3 and is specific of interacting quintessence models. The effective 
equation of state parameter is always positive, so this point never corresponds to an accelerating universe. This 
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point is either unstable or a saddle point depending on the values of P and A. In the limit /3 —>• 0, this point 
merges with the origin. 

• Point F. This point exists only when /l(/3 — A) > —3/2. The energy density of the scalar field and interaction 
components are non-zero, and the energy density of matter vanishes. This point can correspond to an accelerating 
solution when the parameters lie in the range —2 < ^ < 1, and corresponds to either a stable spiral or a saddle 
node. The region in parameter space where this point is stable is indicated by region F in Fig. Thus we can 
conclude there are points in parameter space where this point describes a late time cosmological accelerating 
attractor solution. 


Point 

Existence 

Weff 

Acceleration 

Stability 

O 

VA,d 

w 

No 

Saddle node 

A- 

VA,d 

1 

No 

Unstable node: P > ^, X > — %/6 

Saddle node: otherwise 

A+ 

VA,d 

1 

No 

Unstable node: P < X < x/G 

Saddle node: otherwise 

B 

A= > 3 

w 

No 

Stable node: 3 < A^ < 24/7 & A/3 > 0 

Stable spiral: A^ > 24/7 & A/3 > 0 

Saddle node: XP < 0 

C 

A^ < 6 

A^-3 

3 

A^ < 2 

Stable node: A^ < 3 & A/3 > (A^ — 3) 

Saddle node: /3 < (A^ — 3)/A 

D 

VA,d 

2/3^ 

3 

No 

Unstable node: /3^ > 3/2 & A < (3/3^ + 3)//3 

Saddle node: otherwise 

F 

3 -f 2/3(d - A) > 0 

/3 

x-p 

-2 < A < 1 

Stable spiral: region in Fig. [2] 

Saddle node: otherwise 


TABLE III. Stability of critical points of Model 1 with w = 0 and a = 1. 





FIG. 2. p-X parameter space indicating the region where the critical points B, C and F are stable. 

Analysing the stability of the fixed points, we see that when P > 0 the late time behaviour of the dynamical 
system is qualitatively the same as the case of the canonical scalar field with no interaction term. Although we now 
potentially have up to two extra critical points arising, they are either unstable or saddle points, and at late times 
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the solution will either limit towards the scalar field dominated point C or the scaling solution B. On the other hand 
when /3 < 0 we have new late time behaviour arising. Depending on the parameter values the late time attracting 
solutions are either the scalar field dominated point C, or the point F, a point in which the energy density of the 
matter, scalar field and interaction term is non-zero. This point can represent an accelerating solution for a wide 
range of parameter values. Note that Point F can in principle represent a solution to the cosmic coincidence problem 
since one can attain an accelerating late time attractor where the energy density of dark energy does not dominate 
completely (accelerating scaling solution). 

Now the phase space of the dynamical system is the subset of defined by —oo < a; < oo, y > 0, —oo < z < 
1 — Hence the phase space is non-compact. In order to plot the phase space, we therefore compactify the 

system by introducing the following new variables 

X = arctanx, Y = arctany, Z = arctanz. (99) 

The phase space is now compact, with X,V,Z now lying in the range —7r/2 < X < 7r/2, 0 < V < 7r/2 and 
—7r/2 < Z < arctan(l — tan^ X — tan^T). This is displayed in Fig. [31 One could have chosen other functions to 
compactify, for instance, X = arctanhx, etc. which leads to the same qualitative results. 

We show phase space diagrams with a few trajectories for two distinct cases, when /? > 0 and /3 < 0 in Fig. [3] and 
Fig. m respectively. In Fig. [3] the trajectories start from either the stiff matter states A± and evolve towards point B, 
which serves as the global attractor for this choice of parameters. This dynamics, similar to the case of the canonical 
scalar field with no interaction term, is well suited to describing the late time phenomenology of our universe as it 
has trajectories describing a decelerated to accelerated transition, corresponding to the dominance of dark energy 
over dark matter at late times. However this model suffers from the same issues as the canonical scalar field at early 
times, since the early time attractors are the points A±, which have effective equation of state lUeff = I which is not 
physically viable on the classical level. In Fig.|4|all trajectories again start at either H±, however they now evolve to 
the point F which is the global attractor, which for this choice of parameters is a solution with iCefi = —5/6. Again, 
for suitable choice of the parameters this model can well describe the late time universe phenomenology, however it 
breaks down at early times. 



FIG. 3. Phase space showing trajectories of the dynamical system (1831) - (I85|) when q: = 1,/3 = 1,A = 2 and rc = 0. Point B is 
the global attractor where the universe behaves as though it were completely matter dominated. 

The case a = 3. 

We will now consider the effects on model I of the above dynamical system when we take a different value of a. In 
general when a I, the number of fixed points of the system increases. For simplicity, we will not consider a = 2 
as this gives us critical points at infinity, so we will analyse the case a = 3, again in a matter dominated universe 
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FIG. 4. Phase space showing trajectories of the dynamical system ((8^-(|85l) when O' = 1, /3 = —5, A = 2 and w = 0. Point F is 
the global attractor, and for this choice of parameters it corresponds to a solution with WeS = —5/6. 


with w = 0. Note that values of a different from one correspond to scalar-fluid couplings non-linear in p. Such 
kind of non-linear couplings can only be built using the variational formalism presented in this work, whereas other 
interacting theories can only couple the scalar filed to p linearly, as in the case of Scalar-Tensor theories for example 
[iM3- Depending on the parameter choices, the system now has potentially up to eight fixed points which are 
displayed in Tab. IIVI 


Point 

X 

y 

z 

O 

0 

0 

0 

A± 

±1 

0 

0 

B 

/T1 

V 2 A 

fz 1 

V 2 A 

0 

C 

A 

a/6 

A-# 

0 

D 


0 

1 2/32 

^ 3 

E 

^/6 

0 

0 

3 

~w 

F 

fz 3 

V 2 A-/3 

U-9+2/3(/J-A) 

v^|A-/3| 

A(A-^)-9 


TABLE IV. Critical points of Model I with ui = 0 and q = 3. 

The existence and stability properties of these critical points, along with the values of WeS are displayed in table IVl 
Qualitatively, the properties of the critical points 0,A±,B and C remain the same as the a = 1 case, the only 
modification being the exact regions of ,5 — A parameter space where stability/instability arises are slightly different. 
However the points D, E and F exhibit new behaviour when a = 3. 

• Point D. This point exists for all values of A and /3. It describes a solution where the kinetic energy of 
the scalar field and the energy density of the interaction term is non-zero, but the energy density of matter 
vanishes. Unlike the a = 1 case, this point can now describe an accelerating solution and is also stable for 
certain regions of parameter space. If /3^ > 7/2 this solution describes an accelerating universe and when /3^ > 3 
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and 2/3A < 2/3^ — 9, this point is a stable node. Hence whenever the point is stable, it describes an accelerating 
solution. 

• Point E. This point only exists when /3 > 3. The potential energy of the scalar field vanishes, yet the scalar 
field’s kinetic energy, along with the energy density of the matter and interaction fields, are all non-zero. This 
solution is a scaling solution where the universe behaves as if it were completely matter dominated, yet the 
universe evolves under the influence of all three present fields. This solution is never a late time attractor, it is 
either an unstable node or a saddle point. 

• Point F. This point only exists when 2/3(/3 — A) > 9. As in the a = 1 case, this point represents a solution 
with both scalar field and interaction term components, with the energy density of the matter field vanishing. 
For certain values of parameter space, this solution is a stable spiral and the global attractor, see Fig [5] The 
solution is accelerating only if—2/7 < A//3 < 1. 


Point 

Existence 

WeS 

Acceleration 

Stability 

O 

VA,/3 

w 

No 

Saddle node 

A- 

VA,/3 

1 

No 

Unstable node: /3 > ^ —a/B 

Saddle node: otherwise 

A+ 

VA,/3 

1 

No 

Unstable node: d < 

Saddle node: otherwise 

B 

A^ > 3 

w 

No 

Stable node: 3 < A^ < 24/7 & /3/A > —2 

Stable spiral: A^ > 24/7 j3/\ > —2 

Saddle node: j3/\ < —2 

C 

A^ < 6 

A^-3 

3 

A= < 1 

Stable node: A^ < 3 & A/3 > (A^ — 9) 

Saddle node: A/3 < (A^ — 9) or A^ > 3 

D 

VA,/3 

6-2/3^ 

3 

> 7/2 

Stable node: 13^ > 3 2/3A < (2/?^ — 9) 

Unstable node: < 3/2 & 2/3A > (2/3^ — 9) 

Saddle node: otherwise 

E 

/3^ > 3 

W 

No 

Unstable node: A//3 > —1/2 

Saddle node: A//3 < —1/2 

F 

2/3(d - A) - 9 > 0 

2A+/3 

A-/3 

-2/7 < 1 < 1 

Stable spiral: region VI in Fig. [5] 

Saddle node or unstable node: otherwise 


TABLE V. Stability of critical points of Model I with w = 0 and a = 3. 

The late time behaviour of this model is much more complex than the case a = 1. The qualitative late time 
behaviour of the phase space can be divided into essentially six regions of /3 — A parameter space. Two of these regions 
have two stable critical points, so the late time behaviour is dependent on the initial conditions. The remaining four 
regions have one stable point, corresponding to a global attractor. The different regions are plotted in Fig.[5]displaying 
the corresponding critical points which are stable in that region. In order to plot trajectories in the phase space, we 
again compactify the phase space using the variables (IM)) . We plot some example trajectories for two different choices 
of the parameters in Fig. Inland Fig. [T] The parameter choices for Fig. [S] correspond to the same parameter choices 
used in Fig. [Sj except now with a = 3 instead of a = 1, to enable us to see the effects of changing a. In Fig. [6] the 
qualitative behaviour does not change much from the a = 1 case. Trajectories start at A_, and end at the global 
attractor which is the point B. Many trajectories now pass near the point ZI, which for the choice of parameters 
used has an effective equation of state WsS = 4/3, which is not physically viable. However in Fig. 0 where /3 < 0, 
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Region 

Fixed points stable 

I 

D«feB 

II 

B 

III 

C 

IV 

C&D 

V 

D 

VI 

F 


FIG. 5. /3-A parameter space indicating the different regions where the critical points B, C, D and F are stable. 


all the trajectories start at the point H, and end at the global attractor which is the scaling solution B. This model 
never contains an early time accelerating solution, such a solution is required to describe an inflationary phase of the 
universe. The early time attractors are either A±^ E or D, with D being the only point that can accelerate for certain 
values of parameter space. Unfortunately in the region in which D describes an accelerating solution, the point is a 
saddle, so does not describe an early time inflationary universe. Requiring special initial conditions, for which the 
universe remains in the vicinity of point D for a sufficiently long time, one can however obtain a transient accelerating 
phase, which can still be used to describe early time inflationary stages. 



FIG. 6. Phase space showing trajectories of the dynamical system (1831) - (1851) when a = j3 = 1, \ = 2 and w = 0. 

The analysis of Model I with a = 3 shows the possible interesting applications at cosmological scales of this 
particular coupling between quintessence and dark matter. One can in fact not only obtain late time accelerating 
attractors and scaling solutions, but also transient inflationary epochs useful for early time applications. Moreover 
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FIG. 7. Phase space showing trajectories of the dynamical system (I83l)-(l85l) when a = 3, /3 = —1, A = 2 and w = 0. 

multiple late time attractors can be found for some values of the model parameters. This is a feature which usually 
does not appear in quintessence models of dark energy and can interestingly be analysed with advanced dynamical 
systems methods such as bifurcation theory. The complexity of Model I prevents a detailed study for every value 
of a. However the analyses we considered for the values a = 1 and a = 3 show that these interacting quintessence 
models can successfully be used to describe the dynamics of dark energy at late times. Furthermore the presence of 
accelerating scaling solutions attracting the phase space trajectories at late times indicates a possible solution for the 
cosmic coincidence problem, as usually achieved introducing a coupling between dark energy and dark matter. 


2. Model II 

Now we will analyse the second model of the dynamical system (I83l) - (l85l) where we assume the interaction energy 
density to take the simple form p-mt = lK(l)p, where 7 is a dimensionless parameter. As before we can define an 
effective energy density and pressure, and this time the acceleration equation now reads 

^ = ^[-(1 + w) + {w- l)x'^ + (1 + w)y'^] (100) 

and thus we can define the effective equation of state parameter in this model to be 

Weff = w — {w — \)x^ — (1 + w)y^ (101) 

which we note is independent of the interaction energy density z. The critical points of this model are displayed in 

Tab. ED 

We see that all the critical points are independent of the parameter 7 . We note that the origin of the phase space 
is not a critical point in this model. The existence and stability properties of the points are displayed in Tab. I VIII 
Depending on the value of A, the system has either four or five critical points: 

• Point A± These two points appeared in model I, and correspond to the scalar field kinetic energy dominated 
points with the effective equation of state reducing to that of a stiff fluid iCefi = 1. No acceleration is present 
at this point. These points are either unstable or saddle points depending on whether the absolute value of A is 
less than v^- 
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Point 

X 

y 

Z 

A± 

±1 

0 

0 

B 

[ 3 (l+-w) 

/i" U (i-i-™)(i-™) 

A 2 - 3 (l + m) 

V 2 A 

V 2 A 


C 

A 

a/6 

U-T 

0 

D 

0 

0 

1 


TABLE VI. Critical points of Model II. 


• Point B Unlike the canonical case when no interaction energy is present, this point exists for all values of A. 
This solution is dominated by the energy density of the scalar field and the interaction energy density, and the 
universe behaves as if it is matter dominated at this point. The Jacobian matrix at this point has one zero 
eigenvalue, and two negative eigenvalues, and hence this point is not hyperbolic, meaning its stability properties 
cannot be determined by simply looking at the eigenvalues. To investigate its stability properties, one could use 
centre manifold theory, however here we will just examine this issue numerically. Based on numerical results, 
and by making a comparison with the corresponding point in the non-interacting scalar field model, we postulate 
that this point is stable if 3(1 + w) < A^, and unstable otherwise. 

• Point C This point only exists when A^ < 6. It corresponds to a scalar field dominated universe, with both the 
matter and interaction energy densities vanishing. It is stable if A^ < 3(1 -I- w), and a saddle point otherwise. 
The effective equation of state parameter is A^/3 — 1 and hence the solution describes an accelerating universe 
if A^ < 2. This point is the well known cosmological accelerated expansion driven by a sufficiently flat scalar 
field potential. 

• Point D This point is entirely dominated by the interaction energy z. It exists for all A, and always lies on 
the boundary of the phase space. It is a scaling solution with iCeff = w, so the universe behaves as if it were 
completely matter dominated at this point. This point is non-hyperbolic as one of the eigenvalues of its stability 
matrix is zero, however the matrix always has one positive eigenvalue so the point is unstable. 


Point 

Existence 

Weft 

Eigenvalues 

Stability 

A- 

VA 

1 

3 — 3w, 3 — 3w, 3 -1- y^A 

Unstable node for A > —y/6 

Saddle node: otherwise 

A+ 

VA 

1 

3 — 3w, 3 — Sw, 3 — y^A 

Unstable node: A < a/G 

Saddle node: otherwise 

B 

VA 

w 

0, -1(1 -w) + ^^{1 - w)(24(l -t w)2 - A2(7 + 9w)) 
-|(1 - w) - “ w)(24(l -1- w)2 - A2(7 + 9w)) 

Non-hyperbolic 

Unstable: A^ < 3(1 -I- w) 

C 

A^ < 6 

\2 

— 1 + — 

^,A2-3(1+«;), 

A^ -3(l-tw;) 

Stable node: A^ < 3(1 -I- w) 

Saddle node: 3(1 -I- m) < A^ < 6 

D 

VA 

W 

0, f(w-l),|(u;-tl) 

Unstable 


TABLE VII. Stability of critical points of Model 11. 


In order to plot the trajectories, we again compactify the phase space using the variables (I99|) as previously done 
for model I. We plot the phase space and some trajectories for three different values of A, assuming a matter equation 
of state u> = 0 in Figs. P and fTOl 

We see that this model qualitatively is almost identical to the case of the canonical scalar field when no interaction 
term is present. As seen in Figs. 13 1^ and ITOl all trajectories start at either of the stiff matter points A± and then 
evolve until they reach the global attractor which is either the point B or C, the former corresponding to an effective 
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FIG. 8 . Phase space showing trajectories of the dynamical system when 7 = 1, A = 1.5 and w = 0. 



FIG. 9. Phase space showing trajectories of the dynamical system when 7 = 1, A = 2 and w = 0. 


dust (scaling) solution, while the latter corresponds to an accelerating universe solution if A^ < 2. The main difference 
between this model and the canonical scalar held is the origin is replaced by Point D, corresponding to an interaction 
dominated universe, and the point B now exists for all values of A. The matter dominated solution of the universe is 
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FIG. 10. Phase space showing trajectories of the dynamical system when 7 = 1, A = 4 and ui = 0. 


thus represented by a phase in which the interaction term between dark energy and dark matter dominates. Though 
this phase can present the same dynamics of the ACDM model in the background evolution, differences might arise 
at the perturbation level. The analysis of cosmological perturbations as well as the comparison with observational 
data within the framework of these interacting dark energy models falls outside the scopes of the present work and 
will be left for future studies. 


VI. DISCUSSION AND CONCLUSION 

The main motivation of this work was to study interacting quintessence models which can be derived from a 
variational approach. Many interacting dark energy models considered in the past were constructed by adding suitable 
interaction terms at the level of the field equations. In many cases these terms could not be rigorously motivated 
from an underlining theoretical framework. Here we showed that one can use a rigorous Lagrangian formulation to 
define interacting quintessence models generally depending on one function /(n, s, </>) mixing the dynamical degrees 
of freedom of the matter fluid with the quintessence scalar field. 

This Lagrangian formulation has been used to build first screening mechanisms for Solar System experiments and 
then interesting models of dark energy coupled to dark matter. As shown in Sec. IIVI new mechanisms capable of 
screening the effects of the scalar field at Solar System scales can be obtained. These solutions generalise the well- 
known chameleon mechanism which can now be derived from a variational approach. As an example of the new 
phenomenology that can be obtained in this manner, in Sec. IIVI we considered a cosh-type of interaction proving 
that the scalar field can be highly massive at Solar System distances while being practically massless at cosmological 
scales. 

The cosmological implications of our approach at the background cosmological level have then been studied in 
Sec. El The main difference between the cosmology of our models and the one of previous models is that the matter 
conservation equation does not change in our case. Therefore, matter will always decay according to p oc in 

agreement with the observations. However, the coupling does affect the Hubble constraint equation and therefore the 
total energy density depends on the coupling. 

We proposed two simple interacting models, both of which have a cosmological three dimensional phase space, 
and can be studied efficiently using dynamical systems techniques. The natural phase space of these models is 
unbounded and we introduced new variables which make this phase space compact. The resulting models contain 
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some very interesting phenomenology which we discussed in some detail. In particular for Model I we found late time 
accelerating solutions where the energy density of dark energy either dominates or scales according to the matter 
interacting energy density. This situation is similar to the one generally arising in standard interacting quintessence 
models and can be used to solve or at least alleviate the cosmic coincidence problem. As an example, Fig. ITT] shows 
the behavior of the effective equation of state for one particular trajectory of model I. We note that at late times 
the universe evolves through a matter dominated epoch followed by an epoch of late time accelerated expansion, or 
simply dark energy. This is exactly the behavior observed for our universe, where a matter to dark energy transition 
must occur at late times. Note that for this particular trajectory the early time evolution does not correspond to the 
one we obtain from observations since, as shown by Fig. 1111 a stiff effective equation of state is always attained with 
even excursions to super-stiff region (weff > !)• This is a common feature of quintessence models with exponential 
potential [1^ where the early time attractor is always represented by a scalar field dominated solution with a stiff 
equation of state, and not by a matter dominated solution or even an inflationary epoch as suggested by observations. 
In these situations the model is taken as an effective description valid at late times but not at early times. 


1.0 
0.5 
W/eff 0.0 

- 0.5 

- 1.0 


FIG. 11. The effective equation of state uiefr of model I with q = 3, w = 0, /3 = 2, and A = 1/5. We see a matter dominated 
epoch (corresponding to point E) for some time, followed by the dark energy domiirated late time attractor (point C). 


A 



Log(a) 



Log(a) 


FIG. 12. The effective equation of state Wes of model I with a = 3, w — 1/3, (3 = 2%/3, and A = 2. After a stiff matter period, 
we see an inflationary era (corresponding to point D) for some time, before ending the late time attractor; the scaling solution 
point B. 

Note however that the dynamics of model I (with a = 3) allows also for possible early time inflationary solutions 
described by Point D. This behaviour is neatly showed in Fig. 12 where the effective equation of state for a trajectory 
passing close to Point D has been plotted for w = 1/3 and some particular values of the model parameters. A transient 
early-time inflationary epoch, with Weff = — 1 and a brief excursion into the phantom regime, is dynamically achieved 
after stiff domination, while Point B is the late time attractor describing a radiation dominated scaling solution, 
and not an accelerating solution. Such a situation cannot be used to describe dark energy, but can interestingly be 
adopted to characterize an inflationary phase followed by a radiation dominated era, avoiding in this way the need for 
reheating, a remarkable result by itself. Strong observational constraints are inevitably expected on such a model since 
the presence of the scalar field would be non negligible during radiation domination. Point B being a scaling solution. 
Nevertheless this simple example shows how the scalar-fluid interaction presented in this work can be employed to 
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produce new interesting phenomenology even at early times. 

For a suitable choice of the model parameters, and for a fine tuning of initial conditions, such an inflationary 
solution might even be connected to a late time accelerating attractor through a transient matter dominated epoch. 
This is an interesting possibility since it would allow for a unified description of dark energy and inflation. An in 
depth analysis about the viability of this solution is outside the scope of the present paper and it will be left for future 
studies. Finally another interesting feature obtained within Model I is the possibility of multiple late time attractors 
which are theoretically and mathematically interesting to study. 

Regarding Model II, at the background level it is almost (qualitatively) identical to the uncoupled standard scalar 
field cosmological models. However one expects differences at the perturbation level since in this model the scalar 
field has a non-negligible role during the matter dominated era, this being described by a scaling solution. 

The framework we introduced to study interacting quintessence models is in fact even wider of the one considered 
here. In the present paper we focused on interactions of the form f{n, s, (j)) where no spacetime derivatives appear. 
However, there is in principle no obstacle to consider functions which also depend on those derivatives. Under a 
theoretical perspective, such models are the next logical case to study and will be investigated in the second part of 
this work |38j . One could even consider radically different models where the entropy density s becomes an increasing 
function of time instead of being a constant. Such models, taking into account the second law of thermodynamics, 
would contain a preferred direction of time and provide possible material for quantum gravity speculations. 

In order to investigate the cosmological implications of all these models in more detail, one has to study their 
cosmological perturbations. This would clarify the issues of structure formation and stability, both of which are 
quite challenging to address in general interacting dark energy models. This is because at the level of perturbations 
spacetime derivatives of the coupling enter the perturbed field equations and in general there is no unique procedure 
to promote the background equations to fully covariant equations . Our approach however makes this conceptually 
much easier since we have an underlying theory from which all equations are derived from first principles. In particular 
the fully covariant generalisations of the Einstein and Klein-Gordon field equations can be easily obtained, implying 
that the perturbations equations can be straightforwardly derived and analysed once an interaction has been specified 
at the level of the Lagrangian. 
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Appendix A: Derivation of Eqs. (I39II and (I41II 


In this appendix we will explicitly show the calculations around Eqs. (1391) and m- We start with the first one of 
Eqs. dSni): 


+ ip + p)U>^U''] 

= + p)-ip+ + ip+ p)U^U,V^U'' 


dp 


dp 


= 


dn 


ds 


d(j) 


(Al) 

(A2) 

(A3) 

(A4) 


where in the line (IA2I) we used = 0, obtained from the covariant differentiation of the constraint = —1, 

and in the line (IA3I) we used Eqs. (ITTl) . i.e. = 0 and V^s = 0. 

The second of Eqs. dSS) is longer to verify. Defining p = p + pint and using again Eqs. m and = 0, one 
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finds 


+ C/^C/.)Va + (p + pWU^) 

= h''^VvP + npU^\7\Ufj, 

= h';v,p + n [U^Vxi'pU^) - U^^U^Vxg] 


= KV. in^-f. 


2C/^V[A(/ii7^]) + U^V^iflUx) - 


dn 


dp 




dp 


dn 


9P, 


2u^V[x{pu^^) + - -h’;yuP 

= 2nC/^V[;,(/iC/^])-/i;:|^V,<^. 

In order to show that the first term in the last line above vanishes, we recall Eq. dMl), namely 


pUfj, + + Sd^fj, + (3AO- a — 0 ■ 


Using this we obtain 


2nU^V [A ) = 2U^ [V[AV^](/j + sV[AV^]0 + V [xPa'^ + /3aV[aVo"^] 
= 2U^ (Va/JaV^o^ - V^/JaVao^) 

= 0 , 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 

(All) 


(A12) 


(A13) 

(AM) 

(A15) 


where in the line (IA13I) we used the fact that covariant derivatives commute on any scalar, and in line (|A14I) we 
applied Eqs. ©-(ini). We thus derive 


xT' 


uX 


which corresponds to the second of Eqs. daH). 

To prove Eqs. (HTl) one repeats the calculations (lAlE 




with instead of obtaining simply 


(A16) 


= 0, 

because p does not depend on (j). Analogously one repeats the calculations (IA5|) - (IA11|) which yields 

VVaT"^ = 2nt/^V[A(/iC7M]) ■ 

Then, using again Eq. (jA12l) . and following the similar computations to (IA13I) - (IA15|) one finds 

K.VxT’'^ = -2nU^V[A(/ii„tC/^]), 

implying the validity of Eqs. dUD. 
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